Water-like Anomalies and Breakdown of the Rosenfeld Excess Entropy Scaling 
Relations for the Core-Softened Systems: Dependence on the Trajectory in 

Density- Temperature Plane 
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We show that the existence of the water-like anomalies in kinetic coefficients in the core-softened 
systems depends on the trajectory in p — T plane along which the kinetic coefficients are calculated. 
In particular, it is shown that the diffusion anomaly does exist along the isotherms, but disappears 
along the isochores. We analyze the applicability of the Rosenfeld entropy scaling relations to 
the systems with the core-softened potentials demonstrating the water-like anomalies. It is shown 
that the validity of the of Rosenfeld scaling relation for the diffusion coefficient also depends on 
the trajectory in the p — T plane along which the kinetic coefficients and the excess entropy are 
calculated. In particular, it is valid along isochors, but it breaks down along isotherms. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 



It is well known that some liquids (for example, wa- 
ter, silica, silicon, carbon, and phosphorus) show an 
anomalous behavior [l)-Q:their phase diagrams have re- 
gions where a thermal expansion coefficient is negative 
(density anomaly), self-diffusivity increases upon com- 
pression (diffusion anomaly), and the structural order of 
the system decreases with increasing pressure (structural 
anomaly) 0, 0] • A number of studies demonstrate water- 
like anomalies in fluids that interact through spherically 
symmetric potentials (see, for example, [7J, Il2i . |25 | and 
references therein). 

It was shown 0, Q that thermodynamic and kinetic 
anomalies may be linked through the excess entropy. In 
particular, in Refs. 0, @ the authors propose that the 
entropy scaling relations developed by Rosenfeld can be 
used to describe the diffusivity anomaly. 

In 1977 Rosenfeld proposed the relations connecting 
transport properties of a liquid with the excess entropy 
Q . In order to write down these relations one should use 
the reduced forms of the transport coefficients: 
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where D is the diffusion coefficient. According to the 
Rosenfeld suggestion, the reduced diffusion can be ex- 
pressed in the form 
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where S ex = (S — Sid)/ '(Nks) is excess entropy of the 
liquid and a and b are the constants which depend on 
the studying property Q. The coefficients a and b show 
an extremely weak dependence on the material and can 
be considered as universal. 

In his original works Rosenfeld considered hard 
spheres, soft spheres, Lennard- Jones system and one- 
component plasma [1, 0] . After that the excess entropy 
scaling was applied to many different systems including 




, liquid metals 13l 14 1. 



uids [17j, |18| , network- 
chain fluids [2fJ and 



core-softened liquids 
binary mixtures \l~ 
forming liquids 
bounded potentials [12|, [21 , 

Nevertheless, controversies still remain. For example, 
up to the moment it is not clear whether the Rosenfeld 
scaling relations are applicable to the core softened sys- 
tems. Some publications state that the scaling relations 
are valid for such systems @, [ll| , while in our recent work 
it was shown that the scaling relations may break down 
for the core softened systems [12]. This article presents a 
discussion of this contradiction. Basing on the molecular 
dynamics simulations of the two core-softened systems we 
show that the existence of the water-like anomalies and 
the validity of the Rosenfeld scaling relations depend on 
the trajectory in the (p—T) plane, along which the kinetic 
coefficients are calculated. For the first time it is explic- 
itly shown that the water-like anomalies in the systems 
with core-softened potentials do exists when the kinetic 
coefficients are calculated along the isotherms and do not 
exist along the isochores. Consequently, the exponential 
functional form of the relation between the excess entropy 
and the reduced diffusion coefficient holds along isochors 
while along isotherms one observes its breakdown. 

Two systems are studied in the present work. The first 
one is a repulsive shoulder system (RSS) introduced in 
our previous works 24[ 25|. This system has a potential 



U(r) 



14 



+ n £ ■ I 1 - tanh(/c {> - cti})], (3) 



where a is the "hard"-core diameter, o\ = 1.35 is the 
soft-core diameter, and ko — 10.0. In Ref. [Hj it was 
shown that this system demonstrates anomalous ther- 
modynamic behavior. In our previous publication 12] 
the Rosenfeld relation for this system was studied. It 
was shown that the scaling relation for the diffusion co- 
efficient breaks down for this system in the anomalous 
diffusion region along the isotherms. 
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The second system studied in this work is the core- 
softened system introduced by de Oliveira et al .23 r ] . This 
system is described by the spherically symmetric poten- 
tial represented by a sum of a Lennard- Jones contribution 
and a Gaussian-core interaction (LJG): 
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with a = 5.0, r^ja — 0.7 and c = 1.0. This model can 
qualitatively reproduce water's density, diffusivity, and 
structural anomalies. The diffusivit y o f this system was 
studied, for example, in the papers [Til l23j]. 

In this paper we use the dimensionless quantities: f = 
r/cr, P = Pa 3 /e, V = V/Na 3 = 1/p, f = k B T/e. Since 
we use only these reduced units we omit the tilde marks. 

The simulation setup of the RSS was described in detail 
in Ref. [l2| . The following isotherms are simulated: T = 
0.2; 0.25; 0.3; 0.35; 0.4; 0.5; 0.6; 0.7 and 0.8. 

For the investigation of the LJG potential we simu- 
late a system of 1000 particles in a cubic box for the 
densities ranging from p = 0.01 till p = 0.35 with the 
step Sp = 0.01. The time step used is dt = 0.001. The 
equilibration period consists of 1 • 10 6 time steps and the 
production period - 2.5 • 10 6 time steps. During the equi- 
libration the temperature is kept constant by velocity 
rescaling while during the production cycle NVE-MD is 
used. The equations of motion are integrated by velocity- 
Verlet algorithm. The following isotherms are simulated: 
T = 0.2; 0.3; 0.4; 0.5; 0.6; 1.0 and 1.5. 

The excess entropy in both cases was computed via 
thermodynamic integration method. For doing this we 
calculate excess free energy by integrating the equa- 
tion of state along an isotherm: jjjrpf = jvfcfr — 
_2_ jp p( P )-p k B T ^ rpj^ excesg en t r opy is computed 

vfa <f - U ~ F ^ 
via s ex - NkBT . 

Fig. QJa) presents the diffusion coefficient of the RSS 
system as a function of p and T. Fig. [T^b) demonstrates 
the excess entropy as a function of p and T. From these 
figures it follows that the diffusion coefficient and the 
excess entropy along an isotherm have similar qualitative 
behavior. From Figs, [lja) one can see that there is a 
diffusion anomaly along the isotherms, but there is no 
anomaly along the isochores. Moreover, it is possible to 
find the trajectories in the (p — T) plane, along which 
the diffusion anomaly does not exist. For example, it 
may be shown that there is no diffusion anomaly along 
the isobars. This issue will be discussed in detail in the 
subsequent publication. 

In our previous publication (l~2T | we showed that for this 
system the Rosenfeld scaling relation along isotherms is 
not applicable: the curves demonstrate the self crossing 
loops (see Fig.[2](b)). Fig. [2] (a) represents the logarithm 
of the reduced diffusion coefficient along a set of isochors 
for the RSS. As it can be seen from this figure, the depen- 
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FIG. 1: (Color online) (a) The diffusion coefficient of the RSS 
system as a function of p and T. (b) The excess entropy of 
the RSS system as a function of p and T. 



dence of D* on S ex is linear. However, the slope of the 
line shows an isochor dependence. Fig. [2] (a) shows that 
the slope remains approximately constant for low densi- 
ties (p — 0.3 — 0.55) while on increasing the density the 
slope also increases. From Figs. [2] ((a) and (b)) one can 
conclude that the Rosenfeld scaling relation hold along 
the isochores, but breaks down along the isotherms. 

In order to show that the dependence of the water-like 
anomalies and the Rosenfeld scaling relations on the tra- 
jectory in the (p — T) plane is a general property of the 
core-softened systems we study another core softened po- 
tential introduced above - the LJG potential. Fig. [3Ja) 
presents the diffusion coefficient of the LJG system as a 
function of p and T. As in the case of RSS, one can see 
that at low temperatures there is the diffusion anomaly. 
Fig.[^b) demonstrates the excess entropy as a function of 
p and T. From Figs.[3U(a) and (b)) one can see that there 
is a diffusion anomaly along the isotherms, but there is 
no anomaly along the isochores (compare Figs. [1] ((a) and 
(b)). From these figures it follows that the diffusion co- 
efficient and the excess entropy along the isotherms have 
similar qualitative behavior. However, the location of the 
extremum points of diffusivity and excess entropy along 
the isotherms is different (Fig. [4}. It means that there 
are some regions where one function increases while an- 
other one decreases and vice versa. Clearly, this kind of 
behavior can not be consistent with the Rosenfeld scaling 
formula. From this one can conclude that the Rosenfeld 
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FIG. 3: (Color online), (a) The diffusion coefficient of the 
LJG system as a function of p and T. (b) The excess entropy 
of the LJG system as a function of p and T. 
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FIG. 2: (Color online), (a) Reduced diffusion logarithm for 
RSS along a set of isochors. (b) The logarithm of the re- 
duced diffusion coefficient as a function of S ex along a set of 
isotherms for RSS. 



scaling relation is not applicable to the diffusivity along 
an isotherm for the LJG model. It should be noted that 
the behavior depicted here is consistent with simulation 
results for other water-like fluids (see, for example Ref. 
therefore one can expect that the violation of the Rosen- 
feld scaling relation along the isotherms is the general 
feature for the systems with the water-like anomalies. 

Fig. [5] (a) presents the logarithm of reduced diffusion 
coefficient versus excess entropy along a set of isotherms. 
One can see that at low temperatures the curves demon- 
strate the selfcrossing loops like the ones observed for the 
RSS model (see Ref. [13 anf Fig. [1(b)). These loops be- 
come less pronounced with increasing the temperature 
and at T = 0.5 the self crossing disappears. 

It can be seen from Figs. |3{(a) and (b)), that both 
the diffusivity and the excess entropy along isochors are 
monotonous. It allows to expect that, as in the case of 
RSS, the exponential relation between the diffusion coef- 
ficient and the excess entropy holds along isochors. Fig. [5] 
(b) shows the \n(D*) versus S ex along a set of isochors. 
From this figure one can see that all curves with good ac- 
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FIG. 4: The densities of maximum and minimum of the dif- 
fusivity and the excess entropy for the several isotherms. 



curacy correspond to the Rosenfeld scaling relation. This 
result is in agreement with Refs. Id, 111 |. However, Fig. [5] 
(a) contradicts to the Refs. [H, [Tlf T For example, in Fig. 
3 of Ref. [ll| the authors present the anomalous behavior 
of the diffusion coefficient along the isotherms, but show 
the Rosenfeld relation along the isochores and make the 
wrong conclusion that the Rosenfeld-type excess entropy 
scaling is valid for the LJG system. 

In conclusion, in the present article we carry out a 
molecular dynamics study of the two core-softened sys- 
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FIG. 5: (Color online), (a) The logarithm of the reduced 
diffusion coefficient for the LJG model along the several 
isotherms. The insets correspond to the temperatures T = 0.3 
(upper panel) and T = 0.5 (lower panel) in the larger scale, 
(b) The logarithm of the reduced diffusion coefficient for the 
LJG model along the several isochors. 



terns (RSS and LJG) and show that the existence of the 
water-like anomalies in kinetic coefficients in these sys- 
tems depends on the trajectory in p— T plane along which 
the kinetic coefficients are calculated. In particular, it is 
shown that the diffusion anomaly does exist along the 
isotherms, but disappears along the isochores and iso- 
bars. We analyze the applicability of the Rosenfeld en- 
tropy scaling relations to these systems in the regions 
with the water-like anomalies. It is shown that the va- 
lidity of the of Rosenfeld scaling relation for the diffusion 
coefficient also depends on the trajectory in the p — T 
plane along which the kinetic coefficients and the excess 
entropy arc calculated. In particular, it is valid along 
isochors, but it breaks down along isotherms. 
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